Universal decay law in charged-particle emission and exotic cluster radioactivity 
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A linear universal decay formula is presented starting from the microscopic mechanism of the 
charged-particle emission. It relates the half-lives of monopole radioactive decays with the Q-values 
of the outgoing particles as well as the masses and charges of the nuclei involved in the decay. This 
relation is found to be a generalization of the Geiger-Nuttall law in a radioactivity and explains well 
all known cluster decays. Predictions on the most likely emissions of various clusters are presented. 
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The first striking correlation between the half-lives of 
radioactive decay processes and the Q-values of the emit- 
ted particle was found in a-decay systematics by Geiger 
and Nuttall [l| as, 



logTi/2 = aQ Q 



-1/2 



(1) 



where a and b are constants. However, the Geiger-Nuttall 
law in the form of Eq. (fTJ) has limited prediction power 
since the coefficients a and b change for the decays of 
each isotopic series [2| . It may also change within a sin- 
gle isotopic chain when magic numbers are crossed Q. 
Intensive works have been done trying to generalize the 
Geiger-Nuttall law for a universal description of all de- 
tected a decay events @, i, @, 0, 1, H- Here we present 
a truly universal formula valid for the radioactivity of 
all clusters, including a-particles. This will allow us to 
search for new cluster decay modes and to carry out a 
simple and model-independent study of the decay prop- 
erties of nuclei over the whole nuclear chart. 

We thus observe that the Q-value dependence in 
Eq. |l| is a manifestation of the quantum penetration 
of the a-cluster through the Coulomb barrier. But this 
equation ignores the probability that the a-particle is 
formed on the nuclear surface starting from its four con- 
stituent nucleons moving inside the mother nucleus. This 
is the cause of the limitations of the Geiger-Nuttall law 
mentioned above. In general the decay process, ranging 
from proton to heavier cluster radioactive decays, can 
be described by a two-step mechanism (loj . In the first 
step the formation of the particle and its motion on the 
daughter nuclear surface is established. In macroscopic 
models the clusterization process is described by effec- 
tive quantities adjusted to reproduce as many measured 
half-lives as possible. This procedure has shown to be 
very fruitful, providing a guide to experimental searches. 
In the second step the cluster, with the formation ampli- 
tude and corresponding wave function thus determined, 
is assumed topenetrate the centrifugal and Coulomb bar- 
riers @, i, i, i, 0, EH ■ This second step is well understood 
since the pioneering work of Gamow. It is in fact one of 
the pillars of the probability interpretation of quantum 
mechanics Its great importance in radioactive de- 

cay studies lies in the fact that within a given cluster the 



penetrability process is overwhelmingly dominant. This 
explains the great success of macroscopic models in de- 
scribing radioactive decay. 

We intend here to include the cluster formation prob- 
ability as well as the corresponding penetration through 
the Coulomb barrier. We start from the R-matrix de- 
scription of the cluster decay process [13, [12] • This is the 
basis of all microscopic calculations of cluster decay [l3| . 
The corresponding decay half-life is, 
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where v is the outgoing velocity of the emitted particle 
which carries an angular momentum I. R is a distance 
around the nuclear surface where the wave function de- 
scribing the cluster in the mother nucleus is matched 
with the outgoing cluster-daughter wave function. For 
the distance R we will take the standard value, i.e., 
R = Rq(A 1 J 3 +A 1 J Z ) where Ad and A c are the mass num- 
bers of the daughter and cluster nuclei, respectively. 
is the Coulomb-Hankel function and its arguments are 
standard, i.e., p — /ivR/h and the Coulomb parameter is 
X = 2Z c Zde 2 /hv with /i being the reduced mass and Z c 
and Zd the charge numbers of the cluster and daughter 
nucleus, respectively. The quantity F C {R) is the forma- 
tion amplitude of the decaying cluster at distance R. The 
penetrability is proportional to \H^(x, p)\~ 2 ■ Eq. @ is 
valid for all clusters and for spherical as well as deformed 
cases. The ratio Ni = RF C (R) / H+ (R) , and therefore the 
half- life itself, is independent of the radius R [l3| • 

In microscopic theories the formation amplitude is 
evaluated starting from the single-particle degrees of free- 
dom of the neutrons and protons that eventually become 
the cluster. This is generally a formidable task which 
requires advanced computing facilities as well as suitable 
theoretical schemes to describe the clusterization process. 
It is therefore not surprising that the first calculations of 
absolute decay widths (which require a proper evaluation 
of the formation amplitude) were performed after the ap- 
pearance of the shell model. These calculations had lim- 
ited success due to the small shell model spaces that could 
be included at that time. Only later, with better comput- 
ing facilities, the calculated half-lives started to approach 
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the corresponding experimental values. We will not deal 
with microscopic theories here. For details and references 
on this subject, including an historical background, see 
Ref. [3. 

Our aim is to find few quantities that determine the 
half-life. Expanding in these quantities we hope to be 
able to find, at the lowest order of perturbation, an ex- 
pression of the half-life which is as simple as the Geiger- 
Nuttall law but valid in general, i.e., for all isotopic series 
as well as all type of clusters. With this in mind we notice 
that the Coulomb-Hankel function can be well approxi- 
mated by an analytic formula, which for the I = channel 
reads 

H+( X ,p) « (cot P) 1 ' 2 exp [*(/?- sin/3 cos (3)}, (3) 
where the cluster Q-value is Q c = pv 2 /2 and 

Q C R 
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One sees that cos 2 (3 would be a small quantity if Z c Zd 
is large, i.e., for heavy and superheavy systems. In this 
case one can expand the last term in a power series of 
cos/3. By defining the quantities x' = Z c Zd\J AjQ c and 



pf = ^AZ c Z d (A d /3 + A?*) where A = A d A c /(A d + A c ), 
one gets, after some simple algebra, 
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logT 1/2 = a X ' + &p' + log( 



where 



= e 2 7rV2m/(filnlO) and b = 

— 4e-\/2mi?o/(fidn 10) are constants (m is the nu- 
cleon mass). The first two terms dominate the Coulomb 
penetration and o(3) corresponds to the remaining 
small terms. But still the strong dependence of the 
half-life upon the formation probability in the third 
term of Eq. ([5]) has to be taken into account. It is 
very difficult to make a microscopic calculation of the 
formation amplitude F C (R). But we can extract it from 
the experimental half- lives data by using Eq. @, i.e., 



log\RF c (R)\ = ~log 



^Vo + (x,P)! 2 



(6) 

Taking R = 1.2 fm we evaluated the function 
log I RF C (R) I corresponding to a as well as heavier clus- 
ters. We thus found that the formation probabilities of 
a decays are located in the range log \RF C (R)\ = —1.5 ~ 
—0.75 fm -1 / 2 . The stability of the a decay formation 
amplitude explains the success of the Geiger-Nuttall and 
other empirical laws where formation mechanism is not 
explicitly embedded. However, for all observed cluster 
decays, ranging from a to the heavier 34 Si, the formation 
amplitude changes as much as eight orders of magnitude. 

Yet we found that Eq. ([5]) can still be written as a 
simple linear formula which properly takes into account 
the strong dependence of the formation amplitude upon 
the cluster as well as the mother nuclear structure to a 



first order of approximation. This we have archived by 
exploiting the property that for a given cluster No = 
RF C (R) I Hq (x, p) does not depend upon R. Proceeding 
as above one readily obtains the relation, 



\og\RF c {R)\^\og\R'F c {R')\- 



where R' = R'MT + A c*) 

is a value of the radius that 
differs from R. Since for a given cluster any nuclear struc- 
ture would be carried by the terms RF C (R) and R'F C (R') 
in exactly the same fashion, Eq. ([7]) implies that the for- 
mation amplitude is indeed linearly dependent upon p' . 
Therefore one can write, 
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logTi /2 = ax' + bp' + c. 



(8) 



We emphasize here that the coefficient b in this relation is 
different from that of Eq. ([5]). That is, the terms bp' + c, 
which do not depend upon Q c , have to include the ef- 
fects that induce the clusterization in the mother nucleus. 
Moreover, we found that the term log cot (3/u in Eq. (J5]) 
varies only slightly for all the cases investigated below, 
from a minimum of 0.94 to a maximum of 1.2. The ef- 
fects induced by this variation, as well as the higher order 
terms in Eq. ([5]), are to be taken into account by a proper 
choosing of the constants a, b and c. 

Eq. © holds for all cluster radioactivities. We will call 
this relation the universal decay law (UDL). A straight- 
forward conclusion from the UDL is that log T1/2 depends 
linearly upon x' an d p' ■ This to be valid should include 
the Geiger-Nuttall law as a special case. One sees that 
this is indeed the case since p' remains constant for a 
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Below we will 



given a-decay chain and oc Q 
probe these conclusions, and the approximations leading 
to them. 

We will analyze g.s. to g.s. radioactive decays of even- 
even nuclei. We select 139 a decay events from emit- 
ters with 78 < Z < 108 for which experimental data 
are available. We take the data from the latest compila- 
tions of Refs. [H [TJ] and the lists of Refs. [13, HI- For 
the decay of heavier clusters we have selected 11 mea- 
sured events ranging from 14 C to 34 Si In order to 
perform the calculations one has first to determine the 
values of the constants a, b and c. We carried out an 
extensive search of the best values for these free param- 
eters. Using a fitting procedure for the case of a-decay 
we obtained a = 0.4065, b = -0.4311 and c = -20.7889. 
The quality of the adjustment thus obtained can be seen 
in Fig. [1] where the values of logTi / 2 — bp' as a function 
of x' is shown. The UDL reproduces the available ex- 
perimental half-lives within a factor of about 2.2. This 
compares favorably with modern versions of the Geiger- 
Nuttall law 0]. 

A significant deviation of the UDL in the Figure is the 
nucleus 254 Rf, at x' = 132.16 MeV^ 1 / 2 , for which only 
the lower limit of the half-life is available. This nucleus 
has the value T1/2 > 1.5 ms experimentally [ll|. The 
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TABLE I: Coefficient sets of Eq. ([SJ determined by fitting to 
experimental data in a decay (I), cluster decay (II) and both 
a and cluster decays (III), respectively. The last column is 
given by the Coulomb barrier penetration term of Eq. ([5]) with 
Ro = 1.2 fm. 





1(a) 


II (cluster) 


III(a+cluster) 


IV 


a 


0.4065 


0.3671 


0.3949 


0.4314 


b 


-0.4311 


-0.3296 


-0.3693 


-0.5015 


c 


-20.7889 


-26.2681 


-23.7615 





20 



110 120 130 140 150 160 170 

X' (MeV" 1/2 ) 



FIG. 1: (color online). UDL plots for the a decays of even- 
even nuclei with Z = 78 — 118. The straight line is given as 
a X ' + c. 
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FIG. 2: (color online). Same as Fig. [T] but for the heavier 
cluster decays (left panel) and both a and cluster decays (right 
panel) . 



half- life given by the UDL is T1/2 = 42 ms, correspond- 
ing to a branching ratio of b = 0.055%. A more precise 
measurement of this half-life would be a welcome addi- 
tional test of the UDL. 

We will now analyze cluster decay processes by com- 
paring the predictions of the UDL with the experimental 
data corresponding to the decay of even-even nuclei men- 
tioned above [T^ |. Using the parametrization set II in 
Table U we plotted, as before, the quantity log7i/ 2 — bp' 
as a function of \' ■ As seen in the left part of Fig. [2] 
the agreement between experiment and the UDL is ex- 
cellent. The UDL reproduces the available experimental 
half- lives within a factor of about 4.1. 

Finally we consider all decays together, i.e., a as well 
as heavier clusters. Using the parameter set III of TableU 
we obtained the results shown in the right panel of Fig. [5] 
Again the agreement between the UDL and experiment 



is excellent. 

Using the UDL it is straightforward to evaluate the 
half-lives of all cluster emitters throughout the nuclear 
chart if reliable values of the binding energies are pro- 
vided. This we obtain by using the latest compilation 
of nuclear masses [16j . With the Q-values thus obtained 
we have evaluated the decay half-lives of all isotopes in- 
cluded in that compilation by applying the UDL. We thus 
found that in all cases the experimental values lie between 
the ones calculated by using the parameters of the sets I 
and III in Table HI confirming the prediction power of the 
UDL. We also found that nuclei favoring cluster decays 
are mostly located in the trans-lead region. 

In Table Q] we also give the values of the coefficients 
a and b as provided by Eq. ([5]). It can be seen that 
these values are close to the corresponding fitted values, 
confirming that effects induced by log cot (3/v and higher- 
order terms in Eq. ([5|) are small. 

In summary, we have presented in this paper a simple 
formula that provides with great precision the half-lives 
corresponding to cluster decay. The formula is valid for 
all kind of clusters and for all isotopic series, as expected 
since we derived it from the general description of the 
decay half-life. This formula is of a universal validity and 
therefore we call it universal decay law (UDL). There 
are a few exceptions to this feature, in particular the 
alpha-decay of 254 Rf for which only the lower limit of the 
half-life is available. The UDL predicts that this half-life 
should be T]y 2 = 42 ms. A measurement of this number, 
as well as other cases presented in this paper for heavy 
and superheavy nuclei which may be of interest in present 
experimental facilities, would be most welcome to probe 
the extension of validity of the UDL. This law may also 
help in the ongoing search of new cluster decay modes 
from superheavy nuclei. 
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